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Abstract. We give a brief overview of recent work examining the presence of Q--cIusters in light nuclei within 
the Skyrme-force Hartree-Fock model. Of special significance are investigations into a-chain structures in carbon 
isotopes and "'O. Their stability and possible role in fusion reactions are examined in static and time-dependent 
Hartree-Fock calculations. We find a new type of shape transition in collisions and a centrifugal stabilization of 
the 4a chain state in a limited range of angular momenta. No stabilization is found for the 3a chain. 



1 Introduction 

In light nuclei, cluster structures and large deformation 
may be present even in the mean-field model [l]. For ex- 
ample, ^Be has an a-a structure of "superdeformed" na- 
ture, and the degree of deformation is closely related to 
the orbits of valence neutrons in neutron-rich Be isotopes 
M2I3I . Also there have been many discussions on the nature 
of the Oj state of '^C with well-developed 3a structure, 
possibly a linear chain Il4l? j6l7l . Another example is the 4a 
linear chain band starting around the 4a threshold energy 
region in '^O suggested by Chevallier et al. 1 8 1 through the 
data analysis of the '^C -H a — > **Be + ^Be reaction. This is 
supported by theoretical work|9 |; they analyzed the decay 
widths and discussed that the observed states are possibly 
characterized by well-developed 4a structure. 

Concerning the stability of the linear chain states, in 
conventional models, it has been usually studied through 
the analyses of small vibration around the equilibrium con- 
figuration |10|. However, the bending motion is the most 
important path for the transition to the low-lying states 
Hi]. Therefore, it is necessary to calculate the stability in 
a wide wave function space, which covers not only the lin- 
ear chain configurations but also the lower excited states. 
Utilizing a rather large model space, we have analyzed in 
neutron-rich C isotopes that adding valence neutrons in- 
creases the stability of the linear chain states with respect 
to the bending motion |12 13 1. The basic instability still 
remains, though. 

Even if they are fundamentally unstable, the chain 
structures could still appear as intermediary configurations 
in a collision involving ^Be and/or c-particles, and one 
of the aims of this paper is to explore this possibility us- 
ing time-dependent Hartree-Fock (TDHF) methods. It is 
generally acknowledged that the TDHF theory provides a 
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useful foundation for a fully microscopic many-body the- 
ory of low-energy heavy-ion reactions 1 14 1. Earlier TDHF 
calculations of collisions between light nuclei involving 
cluster-like configurations have been made for the study 
of nuclear molecular resonances [ 15 |. However, due to the 
lack of computational resources these calculations suffered 
from numerical imprecision as well as unphysical symme- 
try assumptions, such as collisions being restricted to ax- 
ial symmetry. Current TDHF calculations are performed 
with high numerical precision and with no symmetry as- 
sumptions as well as using modern Skyrme forces. Re- 
cently, we have shown that when TDHF is combined with 
the density-constraint method |16| dynamical potentials 
and ion-ion interaction barriers can be accurately repro- 
duced II17I18I19L Thus an alternative approach would be 
to investigate the formation and stability of the linear chain 
configuration using the fully microscopic and dynamical 
TDHF theory. 

An important mechanism to stabilize the linear chain 
state could be the rotation of the system; a large moment 
of inertia such as in the linear-chain configuration is fa- 
vored if large angular momentum is given to the nucleus. 
If the angular momentum is too high, however, it leads to 
the fission of the nucleus, and there must exist a region of 
angular momentum where the linear chain configuration 
is stabilized. We focus on this mechanism and discuss the 
angular momentum of the linear-chain configuration based 
on cranked Hartree-Fock. 



2 The 3a chain state in nuclear collisions 

In this work we study TDHF collisions which reproduce 
the linear chain configuration and subsequently decay to 
lower-energy configurations of the system. To our knowl- 
edge such mode changes have never been observed in 
TDHF calculations previously and appear to simulate the 
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suggested astrophysical mechanism for the formation of 
'^C nuclei. 




2500<T<4000 fm/c 
T=5000 fm/c 



Fig. 1. Selected density profiles from TDHF time-evolution of the 
'^He+^^Be head-on collision for initial Be orientation angle /3 = 0° 
using the SLy4 interaction. The initial energy is £;..m. = 2 MeV. 



In our numerical calculations we have chosen a Carte- 
sian box which is 40 fm along the collision axis and 
24 fm in the other two directions. Calculations are 
done in 3-D geometry and using the full Skyrme force 
(SLy4) (201 as described in Ref. \2T\. Using different mod- 
ern parametrizations of the Skyrme force all show the same 
phenomena. We have chosen to study two different col- 
lisions leading to the linear chain configuration, first the 
'^He-i-'^Be system and then the triple collision of three "^He 
nuclei, which may be astrophysically much less probable. 
The Hartree-Fock (HF) state for the '^Be nucleus is axi- 
ally symmetric. In Fig. [1] we show three snapshots from 
the long time evolution of the "^He-n^Be collision at zero 
impact parameter and fic m. = 2 MeV. The top panel of 
Fig.[T]shows the linear chain configuration about which the 
system oscillates for times less than 2, 500 fm/c. In partic- 
ular, it is remarkable that the moving clusters do not equi- 
librate while moving inside the linear chain state but re- 
tain their 2p-2n character, where one observes a complex 
quasiperiodic motion with little damping up to this time. 
Around 2500 fm/c the system starts to bend and acquires a 
somewhat triangular shape as shown in the middle panel of 
Fig- El The system still retains its cluster character with the 
center cluster moving off the reaction plane cut shown in 
the figure, but can be clearly observed in volumetric three- 
dimensional movies of the collision process. The bending 
motion, where the center cluster oscillates somewhat per- 
pendicular to the left and right clusters continues for ap- 
proximately 1000 fm/c, with very little damping. Finally, 
at even longer times the system relaxes into a relatively 
more compact shape (bottom panel of Fig. [T]). Such mode 
changes, where the dynamical energy in the longitudinal 
direction is converted to a transverse mode, while the sys- 
tem retains its cluster structure have never been seen in pre- 
vious TDHF calculations albeit this would not have been 
possible in calculations imposing axial symmetry. Even in 
three-dimensional calculations, for an exactly central col- 
lision, the axial symmetry would be preserved under ideal 
theoretical and numerical conditions. Therefore the mean- 



ing of head-on collision (b - 0) should be interpreted to 
have a small dispersion around this value, which facilitates 
the mode change even for exactly central collisions. 
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Fig. 2. Single-particle parities of the neutron states during the 
collision of the ''He-l-^Be system as a function of time at iic.m. = 
2 MeV. 

We were also successful, for the first time, in creat- 
ing a static Hartree-Fock linear chain state orthogonal to 
the ground state in a fully three-dimensional geometry. 
This was achieved by initializing one of the single-particle 
states to be in the s - d shell with positive parity rather 
than in the p state with negative parity. This results in a 
linear chain state similar to the one shown in the top panel 
of Fig.[T] The fact that this state has two positive parity and 
one negative parity state proves that it is exactly orthogo- 
nal to the '^C ground state which has one positive parity 
and two negative parity states. A similar dependence on 
parity was also studied in cluster model calculations |22|. 
In order to relate this observation to the dynamical mode 
changes discussed above we have used the density con- 
straint method to calculate the potential energy and the 
single-particle parities during the TDHF time-evolution. 
In Fig. |2]we show the neutron single-particle parities as 
a function of collision time. What is striking is that the 
combined system initially has the same parity signature as 
the static linear chain state but at the time of bending one 
of the positive parity states starts to decay towards nega- 
tive parity and this decay continues as the system becomes 
closer to the parity signatures of the ground state. Oscilla- 
tions in the numerically calculated parities stem from the 
fact that these are done during the dynamical evolution of 
the system. The proton single-particle parities are almost 
exactly the same as for the neutrons as anticipated. 

In order to relate the observed mode changes more 
closely to the intrinsic energy of the system we have also 
calculated the potential energy of the system as a function 
of the ion-ion separation distance R IfTTI . For the calcula- 
tion of R we have used the hybrid method described in 1 19 1, 
which relates R to the quadrupole moment of the system 
thus making it possible to have a consistent definition of R 
for large overlaps. The calculated potential energy curves 
are shown in Fig. [3] as a function of three alignments of 
the ^Be nucleus with respect to the collision axis labeled 
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Fig. 3. Potential energy curves for the collision of the '*He+'*Be 
system as a function of R for three initial alignments of the Be 
nucleus and at E^ ,n_ = 2 MeV. 
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Fig. 4. Time development of the potential energy for the head- 
on collision of the ''He+^Be and the ''He+''He+"'He systems for 
Ecm. = 2 MeV. 



as angle y6, and for the entire duration of the collision pro- 
cess. Since for the real system the angular momentum of 
^Be is projected all possible alignments of the ^Be nucleus 
needs to be considered. For all of the alignments the com- 
bined system climbs up a shallow potential barrier height 
of approximately 1 .24 - 1 .44 MeV, the lowest barrier be- 
ing that of the /3 = 0" potential, thus making this alignment 
most probable under threshold conditions. For the poten- 
tial energy curve showing the head-on collision (J3 = 0°, 
black curve) we observe that the system initially relaxes 
to a relatively shallow metastable minimum and oscillates 
about this minimum until approximately T - 2500 fm/c 
at which points it slips down the curve towards the second 
configuration as indicated by three alphas in a triangular 
configuration. After spending some time in this configura- 
tion the system further slips down to even lower energy and 
more compact configuration. The potential energy curves 
corresponding to ^Be initial alignment angles of fi - 45° 
and y6 = 90° (red and blue curves, respectively) undergo 
a different behavior, bypassing the linear chain minimum 
but directly going to the triangular and subsequently to the 
compact configuration, the perpendicular energy collision 
attaining the most compact and lowest energy configura- 
tion. It is interesting to note that all of the potential energy 
curves spend some time in the triangular configuration. 

As an alternate collision leading to the same configura- 
tion we have also studied the collision of three "^He nuclei, 
one at rest at the origin of the collision axis and the other 
two on each side boosted towards the center with 1 MeV 
energy. In Fig.|4]we contrast the time dependence of the 
potential energies for the two different collisions. We ob- 
serve that the three "^He collision process spends consider- 
ably longer time (about 6000 fm/c) undergoing quasiperi- 
odic oscillations with very little damping in the linear chain 
configuration before switching to bending and compact 
modes. 

In order to gauge the stability of the linear chain con- 
figuration we have made systematic studies as a function 
of impact parameter and center-of-mass energy, as well as 
a study using heavier Be isotopes to determine the depen- 
dence on neutron number. In Fig. |5] we show the depen- 



dence of the linear chain survival time on the impact pa- 
rameter for the ^He-i-**Be system at /Sc.m. - 2 MeV and 
P - 0°. We observe that as the impact parameter increases 
the survival time rapidly decreases. This decrease naturally 
happens slower (faster) for lower (higher) energies. We 
have also done a similar study for the time spent in the lin- 
ear chain configuration as a function of the center-of-mass 
energy for the "^He-i-'*Be system for = 0° alignment. We 
decreased the energy in steps of 0.1 MeV to find the lowest 
energy for which we form the linear chain configuration (at 
lower energies the nuclei rebound due to Coulomb repul- 
sion). At this energy of 1.3 MeV the lifetime of the linear 
chain configuration increases to about 2875 fm/c. As the 
energy is increased the lifetime decreases gradually. In or- 
der to study the dependence of the linear chain state on the 
neutron number of the Be isotopes we have repeated all 
of the above calculations using a ''Be nucleus instead. The 
calculations were done by using all the time-odd terms in 
the Skyrme interaction appropriate for an odd- A nucleus. 
While we do find an analogous behavior in this study, the 
lifetime of the linear chain state is approximately 30% less 
than that of the corresponding ^Be system. 
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Fig. 5. Time spent in the linear chain configuration as a function 
of the impact parameter b for the ''He+'^Be system at E^^-^ = 
2 MeV and j8 = 0" alignment. 
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3 Cranked 4a chain states 

To discuss the 4a linear chain state in the rotational frame, 
we perform cranked HF calculations. We self-consistently 
calculate the cranked HF equation, given hy 6{H - coJ) = 
0, where H is the total Hamiltonian, oj is the rotational fre- 
quency, and J is the angular momentum around the y axis. 
We also perform TDHF calculations in order to discuss the 
4a linear chain state in a collision situation. 

We represent the single-particle wave functions on a 
Cartesian grid with a grid spacing of 0.8 fm. The grid 
size is typically 24^ for ground states and 32 x 24^ for su- 
perdeformed states. This accuracy was seen to be sufficient 
to provide converged configurations. The numerical pro- 
cedure is the damped-gradient iteration method |23|, and 
all derivatives are calculated using the Fourier transform 
method. 

We take three different Skyrme forces: Sly6 as a recent 
fit which includes information on isotopic trends and neu- 
tron matter 1.24 J . and SkI3 as well as SkI4 as recent fits 
which map the relativistic isovector structure of the spin- 
orbit force f25\. SkI3 contains a fixed isovector part anal- 
ogous to the relativistic mean-field model, whereas SkI4 is 
adjusted allowing free variation of the isovector spin-orbit 
term. Thus all forces differ somewhat in their actual shell 
structure. Besides the effective mass, the bulk parameters 
(equilibrium energy and density, incompressibility, sym- 
metry energy) are comparable. 
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Fig. 7. Coefficient of the rotational energy, fp- 120, calcuiated 
using the cranking method versus the HF iterations with various 
rotationai frequencies oj. The symbols (b), (c), and (d) correspond 
to the density distribution given in Fig. (6] 
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Fig. 8. Coefficient of the rotationai energy, fp' 120, calcuiated us- 
ing the cranking method versus the rotational frequency oj. The 
lines correspond to the different Skyrme forces as indicated. 



Fig. 9. Angular momentum calculated using the cranking method 
versus the rotational frequency. The lines with solid symbols de- 
note the calculated results for the rigid-body moment of iner- 
tia, while the open symbols denote the results for the cranking 
method; both for the Skyrme forces as indicated. 



We here discuss the stability of the 4a linear chain con- 
figuration in the rotating frame for '^O. To this end, we per- 
form the cranked HF calculations with various rotational 
frequencies, u. For the initial wave function, we use the 
twisted 4a configuration, as shown in Fig.|6](a). Note that 
this initial state violates axial symmetry, which facilitates 
the transition of the initial state to low-lying states includ- 
ing the ground state during the convergence process (this 
was demonstrated for the carbon chain states in II12I13I ). 
We calculate the rigid-body moment of inertia, 0, using 
the total nucleon density at each iteration step. We here 
only consider rotation around the y axis. 

We first investigate the convergence behavior of the 
HF iterations. To check this, we calculate the coefficient 
of the rotational energy, given by tp- 120, at each itera- 
tion step. Figure |7] shows the calculated results with var- 
ious rotational frequencies versus the iterations in the case 
of the SkI3 interaction. The initial state with the twisted 
linear chain configuration is not the true ground state of 
the HF model space and the solution changes into the true 
ground state after some large number of iterations; how- 
ever the situation depends on the value of the rotational 
frequency o). In the figure, we see that the rotational fre- 
quencies OJ - 0.5, 1.0, and 1.5 MeW/h (the dashed, dotted, 
and dot-dashed lines, respectively) lead to the true ground 
state. The corresponding density distribution at the itera- 
tion step of 15000 is plotted in Fig.|6](b). The frequency 
o) = 0.0 MeV/S (the solid line) leads to the quasi-stable one 
(see Fig. |6](c)). At around oj - 2.0 MeW/H, we obtain the 
state (the thin solid line) with the 4a linear chain configu- 
ration, as shown in Fig.|6](d), whereas fission occurs above 
those rotational frequencies (the dot-dot-dashed line). 

We next estimate the range of the rotational frequen- 
cies where the 4a linear-chain configuration can be stabi- 
lized. Figure |8] shows the coefficient of the rotational en- 
ergy, h^/20, versus the rotational frequency oj with vari- 
ous Skyrme interactions. We find stable states for the 4a 
linear chain configuration. For the SkI3 interaction, we ob- 
tain the lower and upper bounds of the rotational frequen- 
cies as 1.8 and 2.2 MeV/ft, between these the 4a linear 
chain configuration is stabilized. The values are 1.9 and 
2.2 MeV/n for the SkI4 interaction and 2.0 and 2.1 MeV/ft 
for the SLy6 interaction, respectively. In these frequency 
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Fig. 6. Total nucleon density distribution calculated using the cranking method for (a) the initial wave function, (b) the ground state, (c) 
the quasi-stable state, and (d) the 4a linear chain state. The isolines correspond to multiples of 0.02 fm/c We normalize the color to the 
density distribution at the maximum of each plot. 



regions where the linear chain configuration is stabilized, 
we can define the rigid-body moments of inertia, which is 
calculated as 0.065 MeV for the SkI3 and SkI4 interactions 
and 0.06 MeV for the SLy6 interaction. There values are in 
very good agreement with the conventional cluster model 
calculations 1261 . 

We also estimate the corresponding angular momen- 
tum where the 4a linear chain configuration is stabilized. 
We calculate the angular momentum using the obtained 
rigid-body moment of inertia and compare those with that 
calculated by the cranking method. The angular momen- 
tum with the rigid-body moment of inertia, Jm, is calcu- 
lated as /rid = Ot^- The angular moment calculated using 
the cranking method, /era, is given by /era J >, where 
< / > is the expectation value of the angular momentum 
in the cranking equation. Figure |9] shows the angular mo- 
mentum thus obtained versus the rotational frequency. In 
the figure, we see that the calculated angular momentum 
using the rigid-body moment of inertia agrees well with 
that of the cranking method, indicating that the rigid-body 
approximation is reasonable for the 4a linear chain states. 
We find that the lower and upper bounds of the angular 
momentum where the 4a linear chain configuration is sta- 
bilized are about 13 and 18 ^ for the SkI3 interaction, 14 
and 18 for the SkI4 interaction, and 16 and 18 for the 
SLy6 interaction, respectively. With such high angular mo- 
menta, very exotic configuration of 4a linear chain can be 
stabilized, however, fission occurs beyond this angular mo- 
mentum region. 

We analyzed such stability of linear chain configura- 
tions of a clusters for three-a in '^C and three-a with va- 
lence neutrons in ^"C. However it is rather difficult to ob- 
tain stable regions of rotational frequency and angular mo- 
mentum in these cases, and stability shown in '^O is con- 
sidered to be a unique character of the 4a configuration. 
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Fig. 11. Coefficient of the rotational energy, fp- 120 calculated 
using the TDHF method for two situations: E = 3.3 MeV with 
b = 8fm, corresponding to / = 2.6S, and E = 18 MeV for the 
same impact parameter with J = I4h. 

shows the starting configuration; the initial velocities are 
in the x-direction. Upon contact the two nuclei are strongly 
attracted and form a compound system that rotates while 
also undergoing strong vibrations as indicated by the other 
snapshots. 

To illustrate the behavior of the moment of inertia we 
show its time development for two typical cases in Fig.fTTI 
The moment of inertia is calculated from the instantaneous 
density as the rigid-body value. They correspond to two 
different angular momenta at the same impact parameter 
of 8 fm, which is selected to get the system to coalesce 
into a strongly deformed rotating configuration; the curves 
also show the quite large oscillations caused by the shape 
changes. Yet the tendency appears to follow that of the 
cranked calculations: for the small angular momentum it 
goes to much smaller moment of inertia, while for J = 14?! 
the system spends quite a long time oscillating around the 
moment of inertia found in the cranking calculations be- 
fore relaxing to a slightly smaller value. 



3.1 Four-a linear chain in collisions of two ^Be with 
TDHF 

To check the importance of this chain configuration as 
an intermediate state in a collision situation, we also per- 
formed time-dependent Hartree-Fock calculations for the 
^Be-i-**Be reaction in a configuration where the axes of the 
collision partners are aligned. Fig. [TOl shows a number of 
typical snapshots of the collision for an impact parameter 
of 6 fm and a relative energy of 2 MeV. The first picture 



4 Summary 

In this work we have performed microscopic dynamical 
calculations of nuclear collisions to study the formation 
of a metastable linear chain state of '^C. The time evolu- 
tion as calculated using the TDHF equations shows a char- 
acteristic quasiperiodic exchange of alpha-like clusters in 
the density function coiTesponding to a quasiperiodic mo- 
tion along a static Hartree-Fock potential energy surface. 
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Fig. 10. Snapshots of the total nucleon density calculated using the TDHF method for the *Be+*Be collision. The impact parameter is 6 
fm and the relative energy is 2 MeV. We normalize the color to the density distribution at the maximum of each plot. 



which is studied using the density constraint procedure. 
We have shown that the calculations lead to the forma- 
tion of a metastable linear chain state of three alpha-like 
clusters which subsequently decays to a lower energy tri- 
angular alpha-like configuration before acquiring a more 
compact shape. This is the first observation of such mode 
changes in TDHF calculations and the results seem to be 
analogous to the suggested astrophysical mechanism for 
the formation of '^C nuclei. 

For the Sa-chain configuration of '^O we found a stabi- 
lization due to centrifugal forces for a range of angular mo- 
menta not including zero. This is a new situation in cranked 
mean-field calculations. An investigation of the presence 
of a chain structure in collisions of ^Be-n^Be shows a ro- 
tating compound system of roughly the right deformation, 
but with superimposed strong shape oscillations. 
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